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Abstract 

We investigate the problem of power flow and its implications to the optimization in power networks. To understand 
how to solve these optimization problems, we look at the injection region of power networks. The injection region 
of a network is the set of all vectors of power injections, one at each bus, that can be achieved while satisfying the 
network and operation constraints. If there are no operation constraints, we show the injection region of a network 
is the set of all injections satisfying the conservation of energy. If the network has a tree topology, we show that the 
injection region with voltage magnitude, line loss constraints, line flow constraints and certain bus power constraints 
has the same set of Pareto optimal points as its convex hull. The set of Pareto-optimal points are of interest since 
these are the the optimal solutions to the minimization of a increasing convex function over the injection region. For 
non-tree networks, we obtain a weaker result by characterize the convex hull of the voltage constraint injection region 
^ ■ for lossless cycles, a lossless cycle with a chord and certain combinations of these networks. The convex hull is of 

C^ , interest since they correspond to optimizing linear functions. 



I. Introduction 

Optimal power flow is a classic problem in power engineering. It is usually given as a static subproblem of the 
security constraint unit commitment problem, in the sense that all the network dynamics such as transients and 
generator behaviors are abstracted away ||T|. The objective of the optimal power flow problem is to minimize the 
cost of power generation in a electrical network while satisfying a set of operation constraints. The cost functions 
are generally taken to be quadratic or linear. This problem has received considerable attention since the late 1960's 
121, and many different algorithms have been developed for it. For a comprehensive review the reader can consult 
Is) and the references within. Despite all the efforts, the optimal power flow problem still remains difficult Q. 

The optimal power flow problem is difficult for two reasons. Firstly, the optimization problem is nonlinear since the 
power injected at each of the buses in the network depends quadratically on the voltages at each bus. Secondly, there 
is typically a large number of different types of constraints. For example, each bus might have voltage magnitude 
and power limits, and each transmission line might have thermal constraints and line flow constraints. Due to 
these two reasons, the optimal power flow problem is a non-convex optimization problem with many constraints, 
and therefore challenging to solve. The traditional approach is to tackle the problem using various heuristics and 
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approximations. One widely used method is to use the so called DC flow approximatioiu where all the lines are 
assumed to be lossless, all voltage magnitude are assumed to be fixed, and all angle differences are assumed to 
be small IS). To contrast with the DC flow approximation, the original optimal power flow problem is sometimes 
called the AC problem. The small angle assumption allows the optimization problem to be linearized and solved. 
However, in practice the lines might be lossy, and the angles might not be small. The small angle assumption comes 
from the fact that connected generators cannot have too much angle difference or they would rotate out of phase 
and become unstable. But in the rapidly emerging area of distributed generation, electricity might be generated by 
renewable sources or from storage devices. If these generators do not have rotating parts, the small angle assumption 
becomes unnecessary. In addition, it is often desirable to make the DC approximation around some AC operating 
point (so called hot start DC); therefore an AC power flow solution is still needed. To solve the full AC problem, 
many global optimization heuristics like genetic algorithms are used, and their effectiveness are generally gauged 
by simulations. But these algorithms do not offer any guarantees about performance and do not offer intuition into 
the structure of the optimization problem. 

A new approach to the traditional optimization methods was taken by the authors in Q . They made the surprising 
empirical observation that in many of the IEEE benchmark networks Q the optimal power flow problem has the 
same optimal value as its convex dual. The main theoretical result is that for a purely resistive network and 
quadratic cost functions with positive coefficients, this convex relaxation is tight. In addition, the result still holds 
if the purely resistive network is perturbed by adding a small reactive part. From this and their observations about 
the IEEE benchmarks, ||6l conjectured that the optimal power flow problem is actually convex for general networks. 
Unfortunately this conjecture is not true since there exist many counter examples lH, Q. A natural question arises: 
if the optimal power flow problem is not convex in general, is it convex for some specific class of networks? 
The results ||6| showed that for 'almost' purely resistive networks the problem is convex, but these networks are 
somewhat unrealistic since practical power networks are mostly reactive instead of resistive. An impetus for this 
paper is to look for some more realistic classes of network for which the optimal power flow problem is convex. 

To find out if there are more realistic networks for which the optimal power flow problem is convex, we focus 
on the feasible injection region of a power network since it allows one to think about power flow in a more abstract 
way and is quite useful in understanding the structure of the problem. The feasible injection region is simply the 
feasibility region of the optimal power flow problem, i.e. the set of all vectors of feasible real power injections 
(both generations and withdraws) at the various buses that satisfy the given network and operation constraints. 
For notational convenience, we drop the word feasible and refer to the region as the injection region. Since the 
optimization problem is solved over the injection region, it is useful to understand the geometry of the region. 

Unfortunately, the injection region is not convex in general. Even though the region is not convex, it still has 

'The DC flow approximation name comes from the fact that historically, engineers used analogue devices that approximated power flow with 
current flow. A perhaps more accurate name would simply be small angle approximation, but the term DC flow approximation is standard and 
we follow it here. 
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some desirable properties for optimization. A subset of the injection region of particular interest is the Pareto frontl^ 
When minimizing an increasing function over a set, the optimal solutions are on the Pareto front. Therefore, even 
though the injection region is not convex, if its Pareto front is convex, the optimization problem is still easy. 

The use of injection region is also useful since it decouples the optimization problem from the physics of power 
flow, thus allowing us to have a higher level view that is often beneficial for other problems in optimization, control 
and pricing in power systems. For example, ifTOl showed there is revenue adequacy in the financial transmission 
rights markets if the injection region has a convex Pareto front. A similar observation is made by ifTTl in the context 
of economic dispatch. This result then can be used if the DC flow assumption is made or if the network is such 
that the AC injection region where the above condition is true. This is similarly the case for many of the recently 
proposed demand response algorithms. 

As a starting point, we look at the injection region of a network with no constraints. In this case, we show 
the injection region is simply the upper half space that satisfies the law of conservation of energyjj Therefore, 
the difficult and interesting part is to quantify how the injection region changes once the operation constraints are 
added. 

There are typically four types of operation constraints in a power network: voltage magnitude, thermal loss in 
transmission lines, line flow limits in a transmission line and bus power limits. If the network is a tree, we show 
that the injection region with voltage magnitude, line loss constraints, line flow constraints and certain bus power 
constraints has a convex Pareto front. Precisely, the condition on the bus power constraints is: each bus is allowed to 
have power upper bounds, but two connected buses cannot both simultaneously have lower bounds. An independent 
work |T3l considered the OPF problem for a tree network, although the authors there used the notion of load 
over-satisfaction and did not consider thermal loss constraints. 

After investigating the tree network, we move onto networks with cycles. Ideally, one would like to state an 
analogous result as in the tree network case. However, we could not yet prove such a strong result. Instead, we 
characterize the convex hull of the voltage magnitude constrained injection region if the network is a cycle with 
lossless links, or a lossless cycle with one chord and certain combinations of these networks. 

The optimal power flow problem has received renewed interest due in part to the emergence of demand response 
and renewable energy in the electricity network. The electricity network is made up of two layers: the transmission 
network and the distribution network. The transmission network consists of high voltage lines that connect big 
generators to cities and towns. The distribution network usually consists of a feeder connected to the transmission 
network, and low voltage lines that connected to end consumers. In addition to the line voltages, the two types 
of networks have different topologies. The transmission network is sparse, but irregular, whereas the distribution 
network is configured to be a tree at any one time of operation. Traditionally, the optimal power flow problem is 
only solved in the transmission network, since the demand in the distribution network is fixed and there is nothing 

^A point in a set is called Pareto-optimal if any coordinate cannot be decreased furtlier without increasing at least one otlier coordinate; the 
Pareto front of a set is simply the set of all Pareto-optimal points. 
^This settles an open question in 1101 , 1121 . 
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to optimize. But with the new 'smart grid' operating paradigm, the demands in the distribution network can be 
variable, and there are more and more renewable and distributed generation in the network as well. In addition, since 
the resistive loss in the distribution network is typically much higher compared to the transmission network, DC 
approximations would perform poorly. Therefore, AC optimal power flow on the distribution network is a relevant 
problem and we show the tree topology of the distribution network simplifies the problem significantly. 

The paper is organized as follows. In Section HIl we establish the notations. Section Hill contains the result about 
the network with no operation constraints. Section |IV] contains results concerning trees. Section [V] contains the 
results about non-tree networks, and Section |VI] concludes the paper. 

II. Model and Notations 

Consider an electric network with n buses. Throughout we assume the network is connected. We work with the 
AC power flow model so in general all variables are complex. We use x to denote vectors, and X to denote matrices. 
Following the convention in power engineering, scalars representing voltage, current and power are denoted with 
capital letters. x0y denote the element- wise product between x and y. Given two real vectors x and y of the same 
dimension, the notation x < y denotes component-wise inequality and x < y denotes component-wise inequality 
with strict inequality in at least one component. We denote Hermitian transpose by (•)^ and complex conjugation 
by conj(-). We write X :>= to mean X positive semidefinite. Given a set ^ G M", convhull(^) denote the convex 
hull of A. 

We write i ~ fc if bus i is connected to k, and i oo k if they are not connected. Let Zik denote the complex 
impedance between bus i and bus k, and yik = — = Qik + jbik denotes the admittance between i and k. Note 
Zik = Zki and yik = Vki- The bus admittance matrix is denoted by Y and defined as 

Y.lr.^ yu if « = fc 

-yik if i ^ k ■ (1) 

if i o^ k 

Y is symmetric, but not Hermitian in general. If Y = Rc(Y), we say the network is purely resistive and if 

Y = j Im(Y), we say the network is lossless. Practical lines are mainly inductive so it is sometimes assumed that 
the network is lossless. Let v = (V^i, F2, . . . , F„) £ C" be the vector of bus voltages and i = (/i, /2, ...,/„) G C" 
be the vector of currents, where /^ is the total current flowing out of bus i to the rest of the network. By Ohm's 
law and Kirchoff's Current Law, i = Yv. The complex power injected at bus i is Si = Pi + jQi = VilJ^ 
where Pi is the real power and Qi is the reactive power Pi positive means bus i is generating real power and 
Pi negative means bus i is consuming real power Let p = (Pi, P2, . . . ,Pn) be the vector of real powers and 
p = Rc(v0conj(i)) = Re(v0(Y-^v^)) = Re(diag(vv^Y-^)) where diag(M) is the vector of diagonal elements 
of a matrix M. The resistive loss on a transmission line between buses i and bus k is given by Lik = \Vi — Vkl'^gik- 
The power flowing from bus i to bus k is denoted Pik, and defined as 

P^k^ReiV^\V^-Vk\*y*k)■ (2) 
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Note L,k = Pik+PM- 

III. Network with No Operation Constraints 

To warm up, let us first consider a network with no operation constraints. Since there are no constraints, the 
injection region is defined as 

7' = {peM" :p = Rc(diag(vv^Y^)). (3) 

In this case, the injection region has a simple characterization. 

Theorem 1. Suppose the power network is connected and have n buses with admittance matrix Y. Let V = {p E 
M" : p = Re(diag(vv^Y^)), V G C"} be the injection region of the network. If the network is not lossless, then 
V is given by 

n 

P = {p e M" : ^ f; > 0} U {0}. (4) 

i=l 

Therefore V is the union of the open upper half space of M" and the origin. If the network is lossless, then V is 
given by 

n 

P = {p e M" : y P. = 0}. (5) 



i=l 



Therefore V is a hyperplane through the origin. 



This result is intuitive pleasing since it says if there are no constraints in the network then the injection region 
is only limited by the law of conservation of energy. Conservation of energy gives the bound ^Y^=\ Pi ^ 0' ^"^ if 
the network is not lossless then J2i=i ^i > 0- Theorem [T] states this is the only constraint on the injection region. 
Clearly the region is convex, thus settling the conjecture in ifTOl . lfT2l . The proof of this Theorem is given in the 
Appendix. 

In practice, there would be some operation constraints on the network. For example, the voltages magnitudes at 
each bus are bounded. Figure [2(a)| shows the injection region of a two bus network with fixed voltage magnitudes. 
The region is an ellipse (without the interior). Even in this simple case, we see that the injection region is no longer 
convex. The next section is devoted to the study of the effect of constraints on the injection regions of tree networks 
and their impUcations to optimization problems. 
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IV. Tree Networks 
A. Pareto Front of Injection Region 

In power networks, we are often interested in solving the following optimization problem 

minimize / (Pi , P2 , • ■ • , ^n ) (6a) 

subject toV^<\Vi\<V^ (6b) 

Lik < kk (6c) 

Ptk < Prk (6d) 

P<P,< P^k (6e) 

p = Re(vv^Y^), (6f) 

where f{Pi,P2, . . . , Pn) is the cost function of the real powers; ( l6b] ). ( l6c] i. ( l6dl i. ( l6el i are the constraints corre- 
sponding to bus voltage, line thermal loss, line power flow and bus power respectively; and (|6B is the physical law 
coupling voltage to power In practice, / is usually an increasing function of the power injections. For example, if 
/(Pi, . . . , P„) = Pi + ■ • ■ + P„, then we are minimizing the loss in the network; or if / is quadratic with positive 
coefficients, then we are minimizing the cost of generation. To better understand the optimization problem in ^, 
we look at the feasible injection region, V, defined as 

P = {p e M" : p = Rc(diag(vv^Y^)), Z, < \VU < V, Wi, 

L.,k < hk Vi - k, P,k < P,k Vi ^k,P^<P,< P,yi}. (7) 

Therefore V is the feasibility region of ^. 

Since we are interested in minimizing increasing functions on V, the relevant geometric objects are the Pareto- 
optimal points of V because only they can be the optimal solutions to ^. The mathematical definition is; 

Definition 1. Let A C M". A point -x. E A is said to be a Pareto-optimal point if there does not exist another point 
X G w4 such that ic < x. Denote the set of Pareto-optimal points of A as 0{A) and is sometimes called the Pareto 
front of A. 

Under many circumstances, the Pareto front of the injection region is convex. By this we mean that the V has the 
same Pareto front as convhull(P). Therefore, (|6]l is a convex optimization problem if / is convex and increasing, 
since we may replace the non-convex region P by a convex region convhull(P) and obtain the same solutions. 
Before stating the general result about 0{V) in Theorem |2l it is instructive to use a two bus example to see what 
are the Pareto-optimal points and the effect of various kinds of constraints on them. 

Consider the two bus example in Figure [T] where y is the line admittance. First consider the case where there 
are only voltage constraints. Suppose that \Vi\ = IV2I = 1 per unit. Then P is an ellipse as shown in Figure [2(a)| 
The bold curve represents the Pareto front and we see it is convex. Next, we consider both voltage constraints 
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Fig. 1. Two bus network. 



and the loss constraint Pi + P2 < I for some I. This is presented by intersecting the elhpse by a half plane as 
in Figure |2(b)| and the bold curve is the resulting Pareto front and it is still convex. Next, consider both voltage 




(a) Voltage constrained. (b) Voltage and loss constrained. 

Fig. 2. Voltage constrained and loss constrained injection regions. 



and bus power constraints. In this case, there are several possibiUties, as represented in Figures [3(a)[ |3(b)| and |3(c)| 
In Figure |3(a)[ both bus have power upper bounds, and the Pareto front is convex. In Figure |3(b)[ Pi has upper 
bound, P2 has both upper and lower bounds, and the Pareto front is convex. In Figure |3(c)[ both buses have lower 
bounds, and we see that the convexity of the Pareto front is not preserved. Note that in the two bus case, the line 




(a) Both buses have power upper (b) P2 have both upper and lower (c) Both are lower bounded. Note 
bounds. bounds. convexity is not preserved in this 

case. 

Fig. 3. Three possible cases of the bus power constrained injection region. 



flow constraints in ( |6dl ) correspond to Figure |3(a)| The intuition gained from the two bus example carries over for 
general trees, and the general statement is given in Theorem |2] 
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Theorem 2. Consider a tree network with n buses. Let the injection region V defined as in (|7]). Suppose if i ^ k 
in the network, then P_.j^ = — oo or Pj. = — oo. Then 0{'P) = C'(convhull(7')). 

The condition on the bus power lower bounds means that if two buses are connected, then not both can have a 
tight bus power lower bound. One case of this being satisfied is if all buses only have upper bounds on its power But 
Theorem|2] states it is not necessary to remove all the lower bounds. For example, consider the tree network in Figure 
m Theorem |2] applies if the following sets of buses have power lower bounds: {1, 4, 5, 6}, {2, 3}, {2, 6}, {3, 4, 5}. 




Fig. 4. A3 level tree with 6 buses. 

Proof: To prove the theorem, first we define an optimization problem in term of the injection region. In this 
optimization problem, we want to write every quantity as a quadratic form of the complex voltages. 

The resistive loss on the transmission line between buses i and k can be written as Lik ~ v^G^^v where G^fc 
is a matrix with the (i, i)th entry and the (fc, fc)th entry being gik, and the (i, fc)th entry and the (fc, i)th entry being 
—gik and all other entries being 0. The power flow from bus i to bus k can be written as Pik = v^Ai^v, where 
Aik is a matrix with (i, i)th entry gik, the (i, fc)th entry j{—gik — jbik), the (k, i)th entry ^{—gik + jbik) and all 
the other entries 0. Similarly, let A^ = i(EiY + Y^E^) where E; is the diagonal matrix with 1 at the {i,i)th 
entry and everywhere else. Then the power injected at bus i is given by Pi = v^A^v. 

Consider the following optimization problem 

n 

J = minimize \J CiPt (8) 

1=1 

subject to K, <\Vi\<V^,\/i 

v^GjfcV < kk yi ^ k 
v"A,kV < P,k Vi - A: 
P, < v^A,v < Pi 
p = Re(diag(vv^Y^)). 
The Ci's can be interpreted as the costs of the power generation and (O is a optimal power flow problem with a 
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linear cost function. To expose the potential non-convexity, we can equivalently write it as 

n 

J = minimize N, CiPi (9) 

subject to VJi < W^i < v\, "ii 

Tr(G,feW) < U^ yi^k 
Tr(A,feW) < P,fc yi^k 
P, < Tr(A,W) < P, 
p == Rc(diag(WY^)) 
W;==0 

rank(W) = 1, (10) 

where W = vv^ and the non-convexity enters as the rank 1 constraint on W. Relaxing this rank 1 constraint and 
eliminating p, we get 

Ji = minimize Tr(MW) (11) 

subject to V^^ < Wu <v\,^i 
W^ 

Tr(G,fcW) < /,fc Vi - fc 
Tr(AifeW) < P,fe Vi - k 
P, < Tr(A,W) < P„ 

where M = i(CY + Y^C) and C = diag(ci, . . . , c„)- Note M is Hermitian. 

Geometrically, the relaxation from d?) to ( fTTT i enlarges the feasible injection region to a convex region given by 

Pw = {p : P = Re(diag(WY^)), 1^2 < w^,^ < y] Vz, 

Tr(G,fcW) < h^ Vz ~ fc,Tr(A,fcW) < P.^ Vz - fc, 

P, < Tr(A,W) < P, Vi, W ^ 0}. (12) 

We want to show that the two regions have the same Pareto front. That is, 0(P) = C'(Pw). Since V^n is convex, 
its Pareto front is easily explored. Note in general Pw ^ convhull(P) and the inclusion can be strict. However, if 
V and P\v have the same Pareto front, then so does convhull(P). The authors in ||6l showed that there is no gap if 
the network is purely resistive and all costs positive. Interpreting this in our language, they showed that the Pareto 
front of the injection region of the resistive network is the same as its convex hull. In contrast, our results are based 
on the topology of the network, and do not need to make assumption about the value of the line admittances. 
The proof of the theorem follows from the following claim. 
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Claim 3. Suppose Ci > for all i. Then the optimal solution to dlll l is unique and has rank 1 if for every connected 
pair of buses [i, k) in the network, one of them do not have tight bus power lower bound. 

This claim is a stronger statement then saying J = Ji, it also states that the optimal solution to the relaxed 
solution is unique. Suppose the claim is true. Then since Pw is convex, we can explore its Pareto front by linear 
functions with positive costs 114]. More precisely, a point pw G "Pw is a Pareto-optimal if and only if it is an 
optimal solution to (fTTT i for some positive costs. From the claim, all the optimal solutions are achieved by a W 
of rank 1, therefore they can be achieved by using a voltage vector v. Therefore if p G Vmv is a Pareto-optimal, 
then p G P. Since Pw 3 P, p is also a Pareto-optimal point of P. So 0{V) ^ 0{V^w)- To show the other 
direction, suppose there exists a point p G 0{V) but not in ©(Pw). Then there is a point p G 0{Vw) such that 
p p. But p G 0{V), contradicting the fact p is a Pareto-optimal point of P. Therefore 0{V) C ©(Pw) and thus 
0(P) = ©(P-w)- It remains to prove the claim. 

We are to show that the optimal solution to (fTTT i. W*, is rank 1. We do this through duality theory. The dual of 
(UBis 

n 

maximize ^(A^y ■ - A^VJ - ^ f^ikhk~ 

i—l ir^k 

n 

i~k i=l 

subject to A + ^ fJ.tkGtk + ^{i^ikA.k + Vki^ki) 

+ ^(a,A,) + M^O, (13) 

where A; and A, are the Lagrangian multiplier associated with the voltage upper and lower bounds and A; = Ai — Aj 
and A = diag(Ai, . . . , A„), ^ik are the Lagrangian multiplier associated with the thermal constraints, Vik and Vki 
are the Lagrangian multipliers associated with the flow constraints, and ct^ and g_^ are the Lagrangian multiplier 
associated with the power upper and lower bounds and <t = ct; — CTj. Note ( fTST l is also the due of ([8]) so the gap 
between J and Ji is called the duality gap. 

Let M = Y^ir^kil^ik^ik) + T.rr^k^V'^kMk + ^kiAki) + Y.l=i{(^i^i) + M. Let W* denote the optimal solution 
of (fTTT i and A* the optimal solution of ( fTsT l, by the complimentary slackness condition lfT4l . 

Tr((A*+M)W*) = 0. (14) 

Since both W* and A* + M are positive semidefinite, (O implies that (A* + M)W* = 0. Therefore W* is in 
the null space of A* + M and rank(A* + M) + rank(W*) < n. So to show rank(W*) = 1 it suffices to show 
rank(A* + M) > n — 1. This is done by considering the topology of the network and thus the structure of M. 

Given a n x n matrix A and a graph G with n nodes, we say that A fits G if for i ^ k, Aik = if and only if 
(i, k) is not an edge in G. The values on the diagonal of A are unconstrained. The next lemma from ifTsl relates 
the topology of a graph and the rank of matrix that fits it. 
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Lemma 4 (Theorem 3.4 in II15I ). Let G be a graph that is a connected tree of n nodes. Suppose A is a n x n 
complex positive semidefinite matrix that fits G. Then rank(A) > n — 1. 

We want to apply this lemma to the matrix A* + M. Since A* is diagonal, only M matters and its entries are 
given by 

-\{{Ci + Cfe + /^jfc + I/jfc + Vkl +0-4+ 0-A;)gifc 

Mife = { j^ j(ci- Cfc + /iife - [iui + Oi - (Tk)bik) if i ^ k ■ (15) 

if i oo k 

Therefore if Mik = if bus i is not connected to bus k. For M to fit the network, Mik needs to be nonzero if i 
is connected to k. 

If i ~ fc, for Mifc to be zero we need 

Ci +Ck+ Pik + Vik + Vkl + Oi + (Tfc = (16) 

c^~C}, + fik - Vki + CT; - cTfc = 0. (17) 

Note pik, Vik and Vki are always nonnegative since they are the Lagrange multipliers associated with upper bounds. 
Suppose the bus power lower bound is not tight for bus i, then cji > 0. Adding (fT6] l with iTl\ gives 2ci + jiik + 
2vik + 2ai = 0, this is not possible since Ci > 0. On the other hand, suppose the bus power lower bound is not tight 
for bus k, then ak > 0. Subtracting ([TtT i from ( fTSI l gives 2ck + Pik + 2i^fci + 2(7^ ~ 0, which is not possible since 
Cfc > 0. Therefore M fits a connected tree. Now apply Lemma|4]to the matrix A* +M gives rank(A* +M) > n— 1, 
therefore rank(W*) < L If the problem is feasible, then rank(W*) = L ■ 

B. Reactive Power and Shunt Capacitances 

Sometimes a bus in the distribution network has a shunt capacitor to ground. This is modeled by adding a 
imaginary diagonal to the bus admittance matrix Y. Let Y be the bus admittance matrix for the network with shunt 
capacitors, then Y = Y + j diag(5io, &201 • • • , ^no) where hm is the shunt capacitance to ground of the ith bus. Let 
B = diag(6io, 620, • ■ • 7 ^no)- Note by adding a shunt capacitance the real power output is not changed since 

p = Rc(diag(vv^Y^)) 

= Re(diag(vv^(Y + B)-^)) 

= Re(diag(vv^Y^)) + Rc(diag(vv^B^)) 

= Re(diag(vv^Y^)), 

where the last equality follows from the fact B is a purely imaginary diagonal matrix. The purpose of adding a 
shunt capacitor is to provide reactive power support. 

Reactive support is needed since some buses may have a reactive power constraint in the form of Q < Qi < Qj, 
without the shunt capacitors, the network may not be able to satisfy this constraint. We show here that a similar 
result as Theorem |2] can be stated if the network have reactive power upper bounds. 
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Corollary 5. Consider a tree network with n buses, each with possibly a shunt capacitor to ground. Assume the 
lines are inductive, that is bik < \/i ^ k. Suppose V^i, V.i, • ■ • , l^n,li„ ore the voltage magnitude constraints, 
{lik,i ^ k} are the line loss constraints, Pi,P_iT--,PmP.n '^''^ ^'^^ ^^^ real power bounds and Qi, ■ ■ ■ ,Qn '^''^ 
the bus reactive power bounds. Let the injection region V defined as 

P = {p e M" : p = Rc(diag(vv^Y^)), F, < \V\, < V, Vi, 
P^k + Pu < kk Vi - fc, P, < P < P.Vz, Q, < 0, Vl}. 



(18) 



Suppose if i ^ k in the network, then Pj = — cxd or P_^. = — cxd. Then 0{'P) = ©(convhul^P)). 



The proof follows exactly the same steps as the proof for Theorem |2] and is omitted. If the network has shunt 
capacitors, it ensures that P will not be empty for a greater set of Q/s than the network without shunt capacitors. 

Note that the above result cannot accommodate lower bound constraints on reactive powers. This again can be 
illustrate by the two bus example. Consider the network in Figure [1] Fix the voltage magnitudes at 1 per unit. 
Figure |5(b)| shows the reactive power injection region (all the reactive powers vectors that can be achieved), and 
the bold segments are the segments that satisfies the reactive power constraint Q < Q2 < Q2- The bold segments 



in Figure [5(a)] shows the corresponding injection region. As we can see, the Pareto front of the injection region is 
not convex. Therefore, in general we cannot extend the result to include reactive power lower bounds. 





(a) Real injection region. 
Fig. 5. Impact of reactive power constraints. 



(b) Reactive injection region. 



C. Simulation Results 

We took the networks from the distribution network database in lfT6l . All the example networks are trees. All 
the nodes have fixed power requirements except the feeder. Therefore this does not exactly fit the conditions in 
Theorem [2] But we found that if the equality demand power requirements are changed to inequalities, and the 
objective function is to minimize the power output from the feeder, all the demand power inequalities are meet 
with equality. This is because increasing the load {Pi becomes more negative) would necessarily increase the loss 
in the system, so if the objective is to minimize the loss, then all the demand power inequalities would be at their 
respective upper bounds. 
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V. Non-tree Networks 

Ideally, one would like to generalize the results for trees to networks with cycles. However, this is difficult. We 
state some partial results in this section, and they will be different than the result stated in Theorem |2] in three 
aspects 

• We focus on lossless networks. 

« Only voltage constraints are considered. 

• Instead of looking at the Pareto front of the injection region, we characterize its convex hull. 

Therefore the results in this section are of a weaker flavor than Theorem|2] since we need to assume that the networks 
are lossless and we only consider voltage constraints. Characterizing the convex hull of the injection region is useful 
if we minimize linear objective functions due to the fact that given a set A, minimizing a linear function over A 
and convhull(^) has the same objective values. In practice, often the cost functions are taken to be linear in the 
power injections. If the convex hull of the injection region can be characterized, then the optimization with linear 
functions can be done efficiently. 

Since only the voltage constraints are considered, the injection region is defined as 

P = {p e M" : p = Re(diag(vv^Y^)), K; < \Vl < V, Vi}- (19) 

We can again define a enlarged convex region Vw as 

V^{pe M" :p = Re(diag(W^Y^)),Z^ < Wu < v\ Vi, 

W ^ 0}. (20) 

We have the following theorem 

Theorem 6. Given a network with n buses represented by its bus admittance matrix Y. Let V and Pw be defined 
as in ( fT9] l and ( |20] l respectively. Then if the network is a lossless cycle or a lossless cycle with one chord, then 
convhull(7-') = T'w- 

The next theorem states that joining the basic types of networks in a certain way preserves the characterization 
result. Given two networks G and iJ, the network K is said to be a 1 -connection of G and H if it is possible 
to decompose K into two components K\ and A'2 such that they have only one node in common and no edges 
between them, where A'l is equal to G and A'2 is equal to H. Note by equal we mean that the admittance matrices 
are identical. In particular, if a line in G or iJ is lossless then its corresponding line in K is also lossless. We say 
K is obtained by 1 -connecting G and H. Figure |7] gives an example of a network obtained by 1 -connecting the 
basic networks in Figure |6] 

Theorem 7. Given a network on n nodes with voltage constraints. Then convhull V = V-w if the network is a 
result of repeatedly 1 -connecting a lossless cycle (with a chord) and a tree. 
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It is simple to check if a network has the topology that satisfies the conditions in Theorem [T] Given a network, 
first decompose it into its one connected parts. This corresponding to identifying the 1 -connecting vertices in the 
network, which can be done in linear time. Then one simply check each of the parts to see if they are a tree, or a 
lossless cycle, or a lossless cycle with a chord. 




(a) (b) (c) 

Fig. 6. Three networks: (a) a tree, (b) a cycle and (c) a cycle with a chord. 




Fig. 7. A graph resulting from 1 -connecting the basic components. 

We suspect that a much stronger statement can be made about lossless networks. In particular, suppose F is a 
lossless network and is one of the two topologies in Theorem |6l then we conjecture 

{p : p = Rediag(vv^Y^),F, < \V\, < V, V^} 
={p : p = Rcdiag(vv^Y^), \VU < F, Vi} 
={p : p = Rediag(vv^Y^), \V\, = V^ Vi}. 

The intuition is that in a lossless network, the power flow is driven by the angle difference between the buses, so 
every injection vector should be achievable by adjusting the phases of the complex voltage vectors while keeping 
the magnitudes constant and there would be no holes in the injection region. 

First we prove Theorem|6] This requires that we prove an analogous result about trees first. Consider the following 
lemma 

Lemma 8. Given a tree network with n buses. Let V and V-w be defined as in ( |19t and ( |20| l respectively. Then 
convhull(7') = Vw- 
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and its relaxation 



Proof: To prove this theorem, it suffices to prove that minimizing linear functions over V and T^^v has the 
same optimal objective value for all coefficients lfT4l . So consider the optimization problem 

n 

J = minimize N^CjPi (21) 

i=l 

subject to y, <\Vi\<V^,\/i 

p = Re(diag(vv^Y^)), 

n 

J I ~ minimize \^ CiPi (22) 

i=l 

subject to Yi < Wu < v\, Vi 
W;>=0 
p = Re(diag(WY^)), 

where the costs can be general (no longer constraint to be positive). We show J = Ji for all c^'s. 
The dual of ( |22] | is 

n 

Ji = maximize ^(Z?A» - y%i) (23) 

subject to A + M ;?= 0. 

From Lemma |4] if the costs are such that M is connected, then the optimal solution to (l22l i is rank 1 and clearly 
J = Ji . If M is disconnected, then M can be written as a block diagonal matrix. Suppose there are K connected 
components of M, then M = diag(Mi, . . . , Mi<-). Since the network is a tree, M^ fits the topology of a tree for 
each i. Then ( |22] | and ( l23T l decomposes into K independent primal-dual subproblems, and we may apply Lemma |4] 
to each of them. Let Wj*, . . . , VFJ denote the optimal solutions to each of the subproblems. By Lemma HI they are 
all rank 1 so we can write W* = v*(v*)^ for each i. An optimal solution W* to the original problem is given 

by W* = v*(v*)'^ where v* = 

Now we prove Theorem |6] The approach is the same as in the proof of Lemma [8] That is, we look at JTH . (l22l) 
and their dual ( |23] |. We say a matrix A is lossless if all the off diagonal terms of A are purely imaginary or 0. We 
prove the following lemma 

Lemma 9. Given a graph on n nodes that is either an odd cycle or a cycle with one chord, if A is lossless, positive 
semidefinite and fits G, then rank(A) > n — 1. 

Theorem|6]can be proved from Lemma|9] Suppose the electrical network is lossless and has the topology of an odd 
cycle or a cycle with one chord. The network being lossless means Y is purely imaginary, and M = i (CY+ Y^C) 
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is also purely imaginary since C is real. Suppose that the costs are such that Mik ^ if {i, k) is connected by a 
line in the network. Since A* is diagonal, the dual matrix A* + M is positive semidefinite, lossless and fits the 
network topology. Apply Lemma |9] shows W* is rank 1 . If the cycle is even, we add a chord between two buses, 
and let the admittance of that chord go to 0. Since all the functions in (l22l) are continuous, the optimal solution of 
the network with a chord approaches the network without the chord as the admittance goes to 0. 

If the costs are such that Mik = even if {i, k) is connected in the network, then M either fits a tree or becomes 
disconnected. If M fits a tree, then apply Lemma |4] If M becomes disconnected, then M can be written as a block 
diagonal matrix. If there are K connected components of M, then M = diag(Mi, . . . , M/^-). Since the network is 
a cycle (with a chord), then M^ is either a tree or a cycle for each i. We can apply Lemma |4] or Lemma |9] to each 
component and obtain an optimal solution W* in the same way as in the tree network case. To finish Theorem |6] 
it remains to proof Lemma |9l 

Proof: Given a graph G, the tree- width of G is a number that intuitively captures how close G is to a tree. For 
example, the tree-width of a tree is 1, and the tree-width of a cycle is 2. The rigorous definition and some methods 
of computing the tree-width the reader may consult IfTTl . A graph of tree- width 2 is also called serial-parallel graph 
or a partial-2-tree. The following lemma collects the known results that we need. 

Lemma 10. If G is a cycle of length n, then the minimum rank of real positive semidefinite matrices fitting G is 
n — 2 i7^i/ . More generally, if the graph has tree-width 2, the minimum rank is n — 2 l\15]l , II 19)1 . 

Given a graph G with n nodes and m edges. We construct a bipartite graph derived from G that we call the 
bipartite expansion of G and denote by B{G). B{G) is a bipartite graph with 2n nodes and 2m edges. Label the 
nodes 1, 2, . . . , n, 1', 2' , . . . ,n' with the bipartition being {1, . . . , n} and {!', . . . , n'}. There is an edge between i 
and k' if and only if i ^ k and {i, k) is an edge in G. If G is an odd cycle then B{G) is also a cycle and if G is a 
cycle with a chord then B{G) has tree-width 2 (a subclass of linear-2-trees in the language of |fT9l ). Two examples 
are given in Figures [8] and |9] If G is an even cycle then B{G) is two disconnected cycles, therefore the assumption 
of odd cycle is needed in the Lemma. 



1 




1'' 



13 3' 



(a) G (b) B(G) 

Fig. 8. (a) shows a 3-cycle and (b) shows its bipartite expansion. 

Given a graph G, suppose A is lossless, positive semidefinite and A fits G. We show that the rank of A cannot 
be lower than n — 1. Suppose A has rank r. Then A can be factored as A = Z^Z for some complex matrix r xn 
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(a) G (b) B{G) 

Fig. 9. (a) shows a 4-cycle with a chord and (b) shows its bipartite expansion. 

matrix Z. Let zi, . . . , z„ G C^ be the columns of Z. They satisfy the graph topology condition 

if i 9^ fc 



zf Zfc = < 



7^ if i - fc 



and the lossless line condition 



Re(zf zfc) = Oif i^k. 
From each complex vector we define two real vectors as 



(24) 



(25) 



Re(z,) 
Im(zi) 


y-i = 


Im(zj) 

-Rc(z,) 



Since z,; G C, then Xi,y,; G M^'". By algebra, Re(zf'^Zfc) = x^x^. ~ y\yk and Ini(zf'^Zfc) = xfy^. In terms of 
x's and y's, (l24l i becomes 

if i oo k 



T 



7^ if i - fc 



and dZSl l becomes 



xf Xfc == yf Yfe = if i ^ fc. 



(26) 



(27) 



Define the matrix B to be the 2r x 2n matrix with columns xi, . . . , x„, yi, . . . , y„. By ( |26] | and dZTb -B fits B{G). 
But if G is an odd cycle or a cycle with one chord, applying Lemma [Tol to B{G) gives rank(_B) > 2n — 2. Thus 
2r > 2n — 2 or r > 77, — 1. ■ 

Now we proceed to the proof Theorem |2l Given a network G, we say the matrix A satisfies G if A fits the 
topology of G and Aik is purely imaginary if the line from bus i to bus k is lossless. We have the following lemma. 

Lemma 11. Given two networks G and H with n and m buses respectively, let K be a network obtained by 
1 -connecting G and H, so K has n + m ~ \ buses. If A is a positive semidefinite matrix that satisfies K, then 
rank(A) > n + m — 2. 

From the basic topologies in Theorem |6l we can apply the Lemma [TT| repeatedly to get Theorem |2l A version of 
Lemma nn just about graphs (without considering lossless lines and such) is known in the graph theory community 
ifTSl . 1201 . We give a proof here to show the additional condition of lossless lines does not change the result. 
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Proof: Let G, H and K be networks given in the statement of the Lemma. Label the buses in K to be 
1, 2, . . . , n — 1, n, ri + 1, 77, + 2, . . . , n + m — 1 where the subnetwork induced by 1, . . . , n — 1, n corresponds to G 
and the subnetwork induced by n, rt+l, n+m— 1 corresponds to H. So bus n is the common bus in the 1-connection. 
Suppose A is a (n + m—\) x (n + m — 1) positive semidefinite matrix that satisfies K and has rank r. Then it is 
possible to factor A as A = Z^Z for some r x (n + rn — 1) matrix Z. Let zi , . . . , z„_)-m_i be the columns of Z. Let 
U be the subspace spanned by zi, . . . , z„_i and V be the subspace spanned by z„-|_i, . . . , z„+„j_i. By construction 
of K, there are no lines between the set of buses {1, . . . , n— 1} and {n+1, . . . ,n+m — 1}. Therefore V is orthogonal 
to U. We may write vector z„ as z„ = u+v+w where u G iY, v e V and w is orthogonal to U and V. Let Zg be the 
matrix with columns zi, . . . , z„_i, u and Z/f be the matrix with columns v, z„_)-i, . . . , •Zn+m-i- Let Ag = Z^Zg. 
Since zf^u = zf^Zn for i = 1, . . . , n — 1, Ag equals the matrix formed by the first n rows and n columns of A. By 
the assumption A satisfies K, so Ag satisfies G. Similarly Z|^Z// satisfies H. By the assumption in the Lemma, 
we have rank(ZG) > ?i — 1 and rank(ZH) > jti — 1, so equivalently dimW > n — 1 and dim V >m — l. Since U 
is orthogonal to V and zi, . . . , z„_|-„j_i spans W + V, rank(A) = dimZ// + dim V > (n — 1) + (m — 1) = n + ni — 2. 

■ 
We cannot hope for a general statement about lossless networks. The Appendix includes a numerical example 
of a complete network with 5 buses, and there is a gap between ( 1211 1 and ( |22] |. 

VI. Conclusion 

We studied the problem of power flow in a network and considered the implication to the optimal power flow 
problem. We focused on the injection region and showed how it can be used to understand the optimal power 
flow problem. For unconstrained networks, we showed the injection region is the entire upper half space. For tree 
networks, we showed that the injection region and its convex hull have the same Pareto front when there is voltage 
magnitude constraints, line loss constraints, line flow constraints, and some subset of bus power constraints. We 
also characterized the convex hull of some commonly occurring networks and a certain combination of them. 
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Appendix 
Proof of Theorem[T] 

The following basic lemma from linear algebra is useful ETI . 
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Lemma 12 (Rank Nullity Theorem). Let A be a n x n real symmetric matrix. Let image(^) and ker(A) denote 
the image and kernel of A, respectively. Then diniiniage(j4) + dinikcr(74) = n and iniage(yl) ® kcr(A) = R", 
where © is the direct sum. 

First consider the case where the network is lossless. Then any feasible injection vector must be on the conser- 
vation of energy plane. We need to show that any point on the plane can be achieved. Since the network is lossless 
Y = j Im(Y) where Im(Y) is a n x n real symmetric matrix and each row of Im(Y) sums to by ([T]). Therefore 
Im(Y) is a generalized graph Laplacian matrix where the admittances can be interpreted as weights on the edges. 
By a standard result in graph theory, dimkcr(Im(Y)) = 1 and kcr(Im(Y)) is spanned by the all one's vector 
1. By Lemma [T2I image(Im(Y)) is the linear subspace in K" orthogonal to 1. Let p° be an injection vector on 
the conservation of energy plane, that is X]"=i ^i = 0- Since l^p*' = 0, there is a unique vector v" such that 
Yv" = p" and 1-^v" = 0. Choose the voltage vector v = (— v" + jl), then 

Re(diag((-v" + jl)(-v" + jl)^Y^)) (28) 

= Rc(diag((v"l^ + l(v")^) Im(Y)) (29) 

+ jdiag((vO(vO)^ + ll^)Im(Y))) 

(a) 
= P , 

where (a) follows from the choice of v" and Im(Y) being symmetric. This finishes the proof for a lossless network. 
Next consider the case where the network is lossy. The proof proceeds in two parts, first we show that the 
conservation of energy boundary X^iLi ^« ~ ^ can be arbitrarily closely from above, and then we show the 
injection region is convex. Since the network is lossy, Rc(Y) is a n x n real positive semidefinite Laplacian matrix. 
By conservation of energy, any power injection vector achieved must satisfy X]"=i ^i > if p ^ 0. Let p*' be a 
vector on the conservation of energy plane. We show there is a voltage vector v that achieves a point arbitrarily 
close to p". Since l-^p° = 0, by Lemma fT2l there is a unique vector v" such that Re(Y)v° = p" and l^v*' = 0. 
Let V = (al + iv°) for some a > and the corresponding injection vector p is 

p = Rc(diag(vv^Y)) (30) 

= Rc(diag((al + -v°)(al + -v°)^(Rc(Y) + jlm(Y))) 
a a 

= diag((aill^ + v"l^ + l(vO)^ + ^v"(v")^) Rc(Y)) 
^^ diag(l(v")^ Rc(Y)) + \ diag(vO(vO)^ Rc(Y)) 
':^^diag(l(pOf) + ijdiag(v"(p°)^) 

= pO + i_diag(vO(pOf), 
where (a) follows from 1 G ker(Re(Y)) and Rc(Y) is symmetrical, (6) follows from the choice of v''. We can 



October 6, 2011 DRAFT 



20 



increase a to make p arbitrarily close to p". For example, if we want ||p — p"||oo < e, then choose 



ln"ll Ilv0|l 
a > 

The next lemma states that V is convex. 

Lemma 13. The injection region V as defined in eqn. is a convex set. 

Theorem [T| follows from Lemma [T3] Since the injection region is convex, and the boundary X]"=i Pi ~ ^ can 
be approached arbitrarily closely from above, it includes the open half upper space. In addition the origin can be 
achieved using the all zeros voltage vector It remains to prove the lemma. 
Proof: For a given network with 7i buses represented by Y, define V-y as 

T'^ = {p e R" : p = Rc(diag(vv^Y^)), ||v||2 < F}, (31) 

where ||v||2 = (X]"=i l^lf)^- ^y approaches the unconstrained injection region as V tends to infinity. V— cannot 
have holes since if p G Vy, then ap G Vy for a G [0, 1]. Therefore to prove the convexity of Vy it suffices to 
prove it has convex boundary. Consider the optimization problem 

n 

J = minimize 2, CiPi (32) 

1=1 

subject to ||v||2 < V 



p = Rc(diag(vv"r^)) 



Relaxing and eliminating p, we get 



Ji = minimize Tr(MW) (33) 

4=1 



subject to y~^ Wii < V 



By changing the costs, we are exploring the boundaries of the two regions with linear functions. We want to show 
that all the point on the boundary of the larger region is in fact in the smaller region. 

First we show that for all M there is an optimal W* for ( [33] ) which is rank I. To solve ( [33] ), expand W in terms 

of its eigenvectors, so W = wiwiwf^ + • • ■ ui„w„w^ where w^ is unit norm and X]"=i Wi <V . Then ( [33] l can 



be written as 



n 

minimize 2, WiW^ AIwi (34) 

i=l 
n 

subject to 2_, ^i ^ ^ 



W = ^(m,w,wf)^0. 



i=l 
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By the well known result about Rayleigh quotients ll22l . to minimize any of the terms wf Mwi, the optimal 
w* = mi, where mi is the eigenvector corresponding to the smallest eigenvector of AI. Therefore the optimal 
solution to ( |33] | is W = J2^=i Wimim^ ~ V mim{^ and is rank 1. 

If mi is not unique, since eigenvector are not continuous in the entries of the matrix, we can perturb Y by an 
arbitrarily small amount to obtain a M that has a unique eigenvector corresponding to the smallest value. Note the 
power vector p is continuous in the entries of Y. From uniqueness of mi and the fact there is no gap between 
(|32] l and ( |33T l, the two regions have the same boundary. Taking V to infinity finishes the proof. ■ 

Appendix 
Additional Result for Trees 

In the proof of Theorem |2] we assume that all the costs are positive. This is actually not needed. Observe that 
as long as Ci 7^ for all i, then M is connected and the optimal solution to (l22l l is unique and rank 1. Define a 
function /(Pi, . . . , Pn) to be strict if for each coordinate i, f is either strictly increasing or strictly decreasing. We 
have 

Corollary 14. Given a tree network with n buses represented by its bus admittance matrix Y. Suppose Vi , V_i , . . . , Vn , V_^ 
are the voltage magnitude constraints. Let the injection region V defined as in ( |19l l. Let f be a strict convex function. 
Then the optimal solution of minimizing f over V and minimizing f over convhul^P) will be achieved at the same 
point in V. 
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Appendix 
Example of Lossless Network with Duality Gap 

Consider a network having the following bus admittance matrix 



Y = j 



and let the costs be C = diag(5,4, 3, 2, 1) and 1/^ = 1 Vi. Let W* be the optimal solution obtained by solving 
(I22I 1. One can check rank(W*) = 2. The optimal value of (|22ll is Ji = Tr(MW*). We need to show that there 
is no rank 1 W in the feasible set such that Ji = Tr(MW). This can be done by using a tool from algebraic 
geometry called Grobner basis ||23]| . It is a computation procedure that allows one to check a feasibility of a set of 
real equality and inequalities. We may use it to check that there does not exist a vector v G C''' such that \V\i < 1 Vi 
and v^Mv= Ji. 

This example is purely numerical since it was generated by a random search. It would be more desirable to 
construct an example with duality gap that offers more intuition. 
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